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A NEW TECHNIQUE FOR FINDING EXACT SOLUTIONS OF 
NONLINEAR TIME-FRACTIONAL WAVE-LIKE EQUATIONS 
WITH VARIABLE COEFFICIENTS 

ALI KHALOUTA AND ABDELOUAHAB KADEM 


Abstract. In this paper, we give the exact solutions in terms of Mittag- 
Leffler functions for nonlinear time-fractional wave-like equations with 
variable coefficients by using a new technique, it is a combination of 
natural transform and variational iteration method called natural varia¬ 
tional iteration method (NVIM). This technique enables us to overcome 
the difficulties arising in identifying the general Lagrange multiplier. The 
numerical results shows that this method is highly effective and is prac¬ 
tically suitable for use in these equations. Three numerical examples are 
given to verify the accuracy and efficiency of the proposed technique. 


1. Introduction 

Fractional partial differential equations are widely used in interpretation and 
modeling of many of realism matters appear in applied mathematics and physics 
including fluid mechanics, electrical circuits, diffusion, damping laws, relaxation 
processes, mathematical biology, and so on. Therefore, to seek the solution 
of fractional partial differential equations is an important aspect of scientific 
research. Several mathematical methods including: Adomian decomposition 
method (ADM) [14], variational iteration method (VIM) [18] homotopy anal¬ 
ysis method (HAM) [15] and homotopy perturbation method (HPM) [4], have 
been developed to obtain the exact and approximate analytic solutions. 

Recently, a new option has appeared, includes the combination of Laplace 
transform, Sumudu transform or natural transform with the previously men¬ 
tioned methods to facilitate and improve the resolution speed of fractional partial 
differential equations. Among wich are: Laplace homotopy analysis method [3], 
Laplace decomposition method [7], Laplace variational iteration method [17], ho¬ 
motopy perturbation Sumudu transform method [16], homotopy analysis Sumudu 
transform method [10], variational iteration Sumudu transform method [1], natu¬ 
ral transform homotopy perturbation method [11], natural decomposition method 
[12], homotopy analysis natural transform method [13]. 

In this paper we will suggest a new technique to the search for exact solutions 
of nonlinear time-fractional wave-like equations with variable coefficients. This 

2010 Mathematics Subject Classification. Primary 35R11, 26A33, Secondary 33E12, 35C05. 

Key words and phrases, nonlinear time-fractional wave-like equations, Liouville-Caputo frac¬ 
tional derivative, Mittag-Leffler functions, natural transform, variational iteration method. 

167 



168 


ALI KHALOUTA AND ABDELOUAHAB KADEM 


technique is a combination of two powerful methods, natural transform and vari¬ 
ational iteration method, called natural variational iteration method (NVIM). 
Consider the following nonlinear time-fractional wave-like equations 


^ ^ Qk-\-m 


( 1 . 1 ) 


*d=i 


* j 


dP 


+ Gii{X, t, v)-^^G 2 i{vxi) + H{X, t, v) + S{X, t), 


i=l 


with the initial conditions 


u(X,0) = ao(X), ut(X,0) = ai(X), (1.2) 

where is the Liouville-Caputo fractional derivative operator of order a, 1 < 

a <2. 

Here X = {xi,X2, ...,Xn), Fiij 1 Gii i,j G {1, 2,..., n} are nonlinear functions of 
X, t and V, F 2 ij, G 2 i i, j G {1, 2,..., n} , are nonlinear functions of derivatives of v 
with respect to Xi and Xj i,j G {1, 2,..., n} , respectively. Also H, S are nonlinear 
functions and k,m,p are integers. 

These types of equations are of considerable significance in various fields of 
applied sciences, mathematical physics, nonlinear hydrodynamics, engineering 
physics, biophysics, human movement sciences, astrophysics and plasma physics. 
These equations describe the evolution of erratic motions of small particles that 
are immersed in fluids, fluctuations of the intensity of laser light, velocity distri¬ 
butions of fluid particles in turbulent flows. 

The plan of this paper is structured as follows. In Section 2, we introduce some 
necessary definitions of the fractional calculus theory and natural transform, two 
theorems are introduced with their proofs related of the natural transform of 
fractional derivatives . In Section 3, we give a brief review of variational iteration 
method (VIM). In section 4, we introduce our results to solve the nonlinear time- 
fractional wave-like equations (1.1) with the initial conditions (1.2) by the natural 
variational iteration method (NVIM). In section 5, we present three numerical 
examples to show the efficiency of this technique and we present our obtained 
results graphically. All numerical results are obtained using Matlab (version 
R2016a). Finally, conclusions are drawn in the last section. 


2. Basic definitions 

In this section, we give some basic notions about fractional calculus, natural 
transform and natural transform of fractional derivatives which are used further 
in this paper. 

2.1. Fractional calculus. We give some basic definitions and important prop¬ 
erties of the fractional calculus theory as the Riemann-Liouville fractional integral 
and Liouville-Caputo fractional derivative (see[9]). 

Definition 2.1. A real function > 0, is considered to be in the space 

/i G M if there exists a real number p > p, so that f{t) = Fh{t), where h{t) 
G C ([0, oo[), and it is said to be in the space if G C)j, n G N. 
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Definition 2.2. The Riemann-Liouville fractional integral operator of order a > 
0 of / G ^ > —1 \s defined as follows 




fo ~ ^ fiOdC t>0,a>0, 

fit), a = 0, 


( 2 . 1 ) 


where r(.) is the well-known gamma function. 


Definition 2.3. The Liouville-Caputo fractional derivative of order a > 0 of 
f £ Cfi,n G N is defined as follows 

jn-af{n)u\^ >0, n_l<a<n, 
a = n. 


D^fit) = 


( 2 . 2 ) 


For the Riemann-Liouville fractional integral and Liouville-Caputo fractional 
derivative, we have the following relation 

n—1 


r D“/(t) = f(t)-Y, (0+) , t > 0. 


fc =0 


kV 


(2.3) 


Definition 2.4. (Partial derivatives of time-fractional order) Let n to be 

the smallest integer that exceeds a. The Liouville-Caputo time-fractional deriv¬ 
ative operator of order a G M"*" is defined as follows 

^ fo [x, Ode, n - 1 < a < n, 


Dfv{x,t)=\ r(n —a) 

v^'^\x, t), a = n. 

Definition 2.5. The Mittag-Leffler function is defined as follows 

OO Yj 

Z 


Ea (Z) = Y, 


rr'n r(na-h 1) 


, a G C, Re{a) > 0. 


n=0 


A further generalization of (2.5) is given in the form 


Ea,l3 (z) — ''Y, 


n=0 


T{na -|- /3) 


,a,f3GC, Re{a) > 0, Re(/3) > 0. 


(2.4) 


(2.5) 


( 2 . 6 ) 


For 0 = 1, Ea (z) reduces to 


2.2. Natural transform. Recently, Fethi Bin Muhammed Belgacem et al. [2] 
introduced a new integral transform, called natural transform, which is applied 
to solve an ordinary and partial differential equations. 


Definition 2.6. The natural transform is defined over the set of functions 
^ = |/W/3M,ri,r2 > 0, |/(t)| < Me^i ,ift G (-1)^ x [0,oo)| , 
by the following integral 


1 r+oo 

[fit)] = R'^ {s,u) = — / e~^ fit)dt, s,u G [0, oo). (2-7) 

Jo 

Some basic properties of the natural transform are given as follows [2]. 
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Property 1: The natural transform is a linear operator. That is, if A and fx 
are non-zero constants, then 

AA+ [\f{t) ± ixg{t)] = AAA+ [fit)] ± [git)]. 

Property 2: If f^^\t) is the derivative of function fit) with respect to 
”t” then its natural transform is given by 

, „n „n—(fc+l) 

jv'+ [/"»(*)] = 

k=0 

Property 3: (Convolution property) Suppose F~^is,u) and are the 

natural transforms of fit) and git), respectively, both defined in the set A. Then 
the natural transform of their convolution is given by 

[if * a) (i)] = uF+is, u)G+is, u), 
where the convolution of two functions is defined by 

if*g)it)=[ fi09it-0d^= [ fit-Odi^df,. 

Jo Jo 

Property 4: 


VNU = 


J\f~^ - -^ = -, n = l,2,... 

Property 5: If a > — 1, then the natural transform of t°‘ is given by 


AA+ [t“] = r (a + 1) 


gO+l ■ 


2.3. Natural transform of fractional derivatives. 


Theorem 2.1. If R^is,u) is the natural transform of the function fit), then the 
natural transform of the Riemann-Liouville fractional integral of fit) of order a, 
is given by 

?/^ 

M+lI“f{t)] = -R+is,u). (2.8) 

Proof. The Riemann-Liouville fractional integral for the function fit) as in (2.1), 

can be expressed as the convolution 

= (2.9) 

r(a) 

Applying the natural transform in the Eq. (2.9) and using properties (3) and 
(5), we have 

AA+[r/(f)] = AA+[^r-i*/(t)] =ii^AA+[r-i]AA+[/(t)] 

, .Cl 1 7/^ 

= u - R^is,u) = — R'^is,u). 

The proof is complete. □ 
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Theorem 2.2. Letn G N* and a > 0 be such thatn — 1 <a<n and R'^{s,u) be 
the natural transform of the function f{t), then the natural transform denoted by 
R+is ,u) of the Liouville-Caputo fractional derivative of f{t) of order a, is given 
by 


f \ „a-{k+l) 


U‘ 


a—k 


D^fit) 


J t=o 


( 2 . 10 ) 


N* |D“/(()1 = = -R+(s,u) - ^ 

Proof. Let 

m = /<”!{(), 

then by the Definition 2.3 of the Liouville-Caputo fractional derivative, we obtain 

t 

D^f{t) = ^ 

r(n -a) J 


t 

1 ( 2 . 11 ) 

-a) J 


r(n- 

0 

Applying the natural transform on both sides of (2.11) and using Eq. (2.8), 


we get 


II 

M* ID"fit)] =M+ [r-“9(l)] = —G+(»,»). 

Also, we have from the property (2) 

AA+b(t)] = AA+[/W(t)], 

n "-“I n-(A:+l) 


( 2 . 12 ) 

(2.13) 


A:=0 


Jt=0 


Hence, (2.12) becomes 

AA+[D“/(t)] = 


,n—a / 


„n-(fc+l) ' 


A:=0 


-i?+(s,u) - 


f \ „a— (fc+1) 


.a—k 


k=0 


D’^f{t)\^^^ = R+{s,u), 


— 1 < n — l<Q!<n. 

The proof is complete. 


□ 


3. Variational iteration method 

Idea of variational iteration method depends on the general Lagrange’s mul¬ 
tiplier method [6]. This method has a main feature, which is the solution of 
a mathematical problem with linearization assumption used as initial approxi¬ 
mation or trial function. This approximation converges rapidly to an accurate 
solution. 
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To illustrate the basic concepts of the VIM, we consider the following nonlinear 
differential equation 

Lv{x,t) + Nv{x,t) = g{x,t), (3-1) 

where L is a linear operator and V is a nonlinear operator, and g{x,t) is an 
inhomogeneous term. 

According to the VIM [5], we can construct a correction functional as follows 

Vn+i{x, t) = Vn{x,t) + [ X{C) [Lv{x, ^) + Nv{x, - g{x, C)] d^, (3.2) 

Jo 

where A(^) is a general Lagrangian multiplier, which can be identified optimally 
via the variational theory and integration by parts. The subscript n denotes the 
nth Qj-fjg];. approximation, Vn is considered as a restricted variation ( i.e. Svn = 0). 

So, we first determine the Lagrange multiplier A(^) that will be identified 
optimally via integration by parts. The successive approximations Vn+i,n > 0 
of the solution v{x, t) will be readily obtained upon using the obtained Lagrange 
multiplier and by using any selective function vq . Consequently, the Solution 


v{x, t) 


lim Vn{x,t). 


4. Solution of nonlinear time-fractional wave-like equations by 
the natural variational iteration method (NVIM) 

Theorem 4.1. Consider the nonlinear time-fractional wave-like equations (1.1) 
with the initial eonditions (1.2). 

Then, by the NVIM the exact solution of Eqs. (1.1) and (1.2) is given as a 
limit of the suecessive approximations {vn{X,t)}^yQ , in other words, 

v{X,t) = lim Vn{X,t). 

n —^oo 


Proof. In order to achieve our goal, we consider the following nonlinear time- 
fractional wave-like equations (I.l) with the initial conditions (1.2). 

First we define 


Nv 

Mv 

Kv 


^ Flij {X, t, v) ’ '^^3 ) ’ 


* 3 


dP 


Gu{X, t, v)-^G2i{vxi), 


i=l 


H{X,t,v). 


(4.1) 


Eq. (1.1) is written in the form 


Dfv{X,t) = Nv{X,f) + Mv{X,t) + Kv{X,t) + S{X,t). (4.2) 


Applying the natural transform on both sides of (4.2) subject to the initial 
conditions (1.2) and using the Theorem 2.2, we get 

1 7 / 7 /^ 

= -v{X, 0) + -^vt{X, 0) + —AA+ [5(X, t)] 
uY 

+ —AA+ [Nv{X, t) + Mv{X, t) + Kv{X, t)] . 


[u(X,t)] 


(4.3) 
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After that, let us take the inverse natural transform on both sides of (4.3) we 
have 

v{X, t) = L{X, t) + [Nv{X, t) + Mv{X, t) + Kv{X, t)]^ , (4.4) 

where L(X, t) is a term arising from the source term and the prescribed initial 
conditions. 

Take the first partial derivative with respect to of Eq. (4.4) to obtain 

0 = -v{X,t)--L{X,t) (4.5) 

[NviX,t)+MviX,t)+Kv{X,t)]^ . 

According to the variational iteration method, we can construct a correct func¬ 
tional as follows 


Vn+l{X,t) = Vn{X,t) (4.6) 

t 

-I 

0 

Or 

Vn+i{X,t) = L{X,t) (4.7) 

+AA-1 [Nvn{X,t) + MvniX,t) + KVniX,t)]] . 


- ^ 


dt 


Therefore, the exact solution of Eqs. (1.1) and (1.2) is given as a limit of the 
successive approximations {vn{X,t)}^yQ , in other words. 


v{X,t) = lim Vn{X,t). 


This completes the proof. 


□ 


5. Illustrative Examples and Numerical Results 

In this section, we solve three numerical examples of nonlinear time-fractional 
wave-like equations with variable coefficients by using the NVIM. 

Example 5.1. Consider the 2-dimensional nonlinear time-fractional wave-like 
equation with variable coefficients 

= -T^ivxxVyy) - -v,l<a<2, (5.1) 

with the initial conditions 

v(x, y, 0) = vt{x, y, 0) = (5.2) 

where Df is the Liouville-Caputo fractional derivative operator of order a, v is 
a function of (x, y, t) x M+. 

According to (4.7), we can construct the iteration formula as follows 
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Vn+i{x,y,t) = e^y + te^y+M-^ 

92 


-AA+ 


52 


dxdy 


{XnxxX' 


nyy) 


(5.3) 


dxdy 


{xyVfixVny') ^r? 


Using the iteration formula (5.3), we obtain 

^a+l 


Vo = 

{l + t)e^ 

Vl = 

+ 

1 

V2 = 

Tl -|- t — 




r(a + 1) r(a + 2) 






t 




j.2a 


+ ■; 


+ ■ 


t 


2a+l 


r(a + l) r(a + 2) r(2a + l) r(2a + 2) 




Then, the general term in successive approximation is given by 


Vnix,y,t) = 


n 


k=0 


(_l)fc^fca ^ 


c^y 


^r(A:a + l) r(A:a + 2 ) 

So, the exact solution of Eqs. (5.1) and (5.2) in a closed form is given by 


v{x,y,t) = lim Vnix,y,t) = 


00 


k=0 


^_ -^'^k^ka ^ 


^T{ka + 1) 
= (E„(-r) + fE„,2(-t“))e"T 


+ 


r(A:a + 2) ) 




where £’„(—1“) and Ea^ 2 {—i^) are the Mittag-Leffler functions, defined by Eqs. 
(2.5) and (2.6). 

Eor a = 2, then 

v{x^ y, t) = {E2{—t^) + tE2,2{—t^)) e^y = {cost + sint)e^^. 

which is the same result as those obtained by the NIM and NHPM [8] for the same 
example. The surface and behavior of the solution for this example is graphically 
presented in Figures 1 and 2 for different fractional orders of a. 


(a) (b) 



0 0 0 0 

t X t 


Figure 1. The surface graph of the 3*^ order approximate solu¬ 
tion by NVIM and the exact solution of Example 5.1 when y = 0.5: 
(a) V when a = 1.5, (b) v when a = 1.75, (c) v when a = 2, and 
(d) V is exact. 
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Figure 2. The behavior of the exact solution and the 3*^ order 
approximate solution by NVIM of v for different values of a for 
Example 5.1 when x = y = 0.5. 


t 

a = 1.7 

a = 1.8 

a = 1.95 

a = 2 

exact solution 

I'^exact '^NVI m\ 

0.1 

1.3953 

1.3999 

1.4046 

1.4058 

1.4058 

3.2196 X 10"^^ 

0.3 

1.5522 

1.5735 

1.5991 

1.6061 

1.6061 

2.1569 X 10-9 

0.5 

1.6359 

1.6755 

1.7272 

1.7424 

1.7424 

1.3095 X 10-^ 

0.7 

1.6540 

1.7088 

1.7854 

1.8093 

1.8093 

1.9680 X 10-6 

0.9 

1.6137 

1.6775 

1.7728 

1.8040 

1.8040 

1.4947 X 10-6 


Table 1. The numerical values of the 4*^ order approximate so¬ 
lution and the exact solution for Example 5.1 when x = y = 0.5. 


Example 5.2. Consider the following nonlinear time-fractional wave-like equa¬ 
tion with variable coefficients 

^2 q2 

Dfv = v‘^-^{vxVxxVxxx) + vl^ivlx) “ ISu® + v, 1 < a < 2, (5.4) 

with the initial conditions 

i;(x,0) = e^, vt{x,0) = e^, (5.5) 

where is the Liouville-Caputo fractional derivative operator of order a, and 
u is a function of {x, t) G ]0,1[ x M"''. 

According to (4.7), we can construct the iteration formula as follows 


Vn+i{x,t) = e^ + te^+J\f ^ —AA"*" 


Xji Q^2 i'^XixXnxx'Vnxxx') (5'6) 


^2 

~^^nx g^2 ^^rixx) 4 " 
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Using the iteration formula (5.6), we obtain 

Vo = (l + t)e^, 

( t°‘ \ 

_ / ^a+l ^2a ^2a+l 

“ 1^1 + ^ + r(a + l) + r(a + 2) ^ f{2a+\) ^ r(2a + 2) 


and so on. 

Then, the general term in successive approximation is given by 

^ / ^ka ^ka-\-l \ 

So, the exact solution of Eqs. (5.4) and (5.5) in a closed form is given by 

/ -j-ka ^/cQ+1 \ 

v(x.t) = + 

= (E„(r) + fE„,2(t“))e", 

where Ea{t°‘) and Ea^ 2 {t°‘) are the Mittag-Leffler functions, defined by Eqs. (2.5) 
and (2.6). 

Eor 0 = 2, then 

v{x,t) = {E2{f) + tE2,2{t^)) e" = e^+*. 

which is the same result as those obtained by the NIM and NHPM [8] for the same 
example. The surface and behavior of the solution for this example is graphically 
presented in Figures 3 and 4 for different fractional orders of a. 


W (b) 



Figure 3. The surface graph of the 3*^ order approximate solu¬ 
tion by NVIM and the exact solution of Example 5.2: (a) u when 
a = 1.5, (b) V when a = 1.75, (c) v when a = 2, and (d) v is 
exact. 
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t 


Figure 4. The behavior of the exact solution and the 3*^ order 
approximate solution by NVIM of v for different values of a for 
Example 5.2 when x = 0.5. 


t 

a = 1.7 

a = 1.8 

a = 1.95 

a = 2 

exact solution 

I'^exact '^NVI m\ 

0.1 

1.8357 

1.8298 

1.8236 

1.8221 

1.8221 

4.1350 X 10“^^ 

0.3 

2.2994 

2.2697 

2.2350 

2.2255 

2.2255 

2.7750 X 10-9 

0.5 

2.8800 

2.8174 

2.7402 

2.7183 

2.7183 

1.6907 X 10-^ 

0.7 

3.5940 

3.4901 

3.3585 

3.3201 

3.3201 

2.5543 X 10-6 

0.9 

4.4670 

4.3129 

4.1140 

4.0552 

4.0552 

1.9535 X 10-6 


Table 2. The numerical values of the 4*^ order approximate so¬ 
lution and the exact solution for Example 5.2 when x = 0.5. 


Example 5.3. Consider the following one-dimensional nonlinear time-fractional 
wave-like equation with variable coefficients 

d 

Dfv = x^ — {vxVxx) - X^ivxxf -V, 1 <a<2, (5.7) 

with the initial conditions 

u(x,0) = 0, vt{x,0) = x'^, (5-8) 

where Df is the Liouville-Caputo fractional derivative operator of order a, and 
u is a function of (x, t) (x, t) G ]0,1[ x M+. 

According to (4.7), we can construct the iteration formula as follows 


Vn+l{x,t) =tx‘^ +J\f 


-1 


U 


-AA+ 


X Q^{Vnx'Vnxx') X {Vnxx') ^r? 


(5.9) 


Using the iteration formula (5.9), we obtain 


ta+i 


Vq = tx ,Vi = t — 


r(Q! -|- 2) 


X ,V2 = [t- 


+ l ^20 + 1 

-b 


r(a-b2) r(2a-b2) 


and so on. 
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Then, the general term in successive approximation is given by 

(_l)fc^fca+l\ 


Vn{x,t) = 


\k=0 


T{ka + 2) 


So, the exact solution of Eqs. (5.7) and (5.8) in a closed form is given by 


v{x,t) = lim Vn{x,t) = > 

77 .—\ ^ ^ 


{-Xft 




. 2=0 


T{ka + 2) 


= x^ {tEa,2{-t°‘)), 


where is the Mittag-Leffler function, defined by Eq. (2.6). 

Eor a = 2, then 

v{x, t) = x^ (t£' 2 , 2 (—2^)) = x^ sint. 

which is the same result as those obtained by the NIM and NHPM [8] for the same 
example. The surface and behavior of the solution for this example is graphically 
presented in Figures 5 and 6 for different fractional orders of a. 


(a) (b) 



0 0 0 0 ’ 
X t X t 


Figure 5. The surface graph of the 3*^ order approximate solu¬ 
tion by NVIM and the exact solution of Example 5.3: (a) v when 
a = 1.5, (b) V when a = 1.75, (c) v when a = 2, and (d) v is 
exact. 



Figure 6. The behavior of the exact solution and the 3*^ order 
approximate solution by NVIM of v for different values of a for 
Example 5.3 when x = 0.5. 
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t 

a = 1.7 

a = 1.8 

a = 1.95 

a = 2 

exact solution 

\v exact — VnVIm\ 

0.1 

0.02488 

0.02492 

0.02495 

0.02496 

0.02496 

6.8887 X lO-^*^ 

0.3 

0.07271 

0.07319 

0.07374 

0.07388 

0.07388 

1.3549 X 10“^^ 

0.5 

0.11604 

0.11752 

0.11934 

0.11986 

0.11986 

1.3425 X 10"® 

0.7 

0.15325 

0.15615 

0.15994 

0.16105 

0.16105 

2.7677 X 10"® 

0.9 

0.18327 

0.18777 

0.19394 

0.19583 

0.19583 

2.6495 X 10“'^ 


Table 3. The numerical values of the 4*^ approximate solution 
and the exact solution for Example 5.3 when x = 0.5. 


6. Conclution 

In this paper, a new technique called natural variational iteration method 
(NVIM) is presented for hnding the exact solutions of the nonlinear time-fractional 
wave-like equations with variable coefficients. This technique is used to overcome 
the complexity of identifying the general Lagrange multiplier. When compared 
with the existing published techniques, it is easy to notice that the new technique 
has many advantages. It is straightforward, easy to understand, and fast, requir¬ 
ing much less computations to perform a limited number of steps of the simple 
procedure that can be applied to find the exact solution of a wide range of types 
of fractional partial differential equations. Furthermore, there is no need for using 
linearization or restrictive assumptions when employing this new method. The 
method was tested on three numerical examples on different situations. Numeri¬ 
cal results obtained by the technique confirm the easily, accurately and efficiency 
of the technique, and therefore, can be widely applied in others fractional partial 
differential equations. 
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